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Abstract. We study the spherically symmetric collapse of a perfect fluid using area- 
radial coordinates. We show that analytic mass functions describe a static regular 
centre in these coordinates. In this central singularity can not be realized 

without an infinite discontinuity in the central density. We construct mass functions 
involving fluid dynamics at the centre and investigate the relationship between those 
and the nature of the singularities. 
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1. Introduction 

Gravitational collapse is one of the main issues of gravitational physics. The nature 
of spacetime singularities formed in gravitational collapse has been studied for a long 
time. In these studies, it is commonly assumed some sort of spacetime symmetry, such 
as spherical symmetry. 

The final fate of spherically symmetric dust collapse has been extensively studied. 
In this general exact solution to Einstein's field equations is available. It was 

proved that a collapsing dust ball results in naked or covered singularities pQ. The 
condition for both outcomes has been derived in terms of the initial data [21 EH E] 
and the structure of naked singularities formed in this system has also been studied in 
detail jSE]. 

As for more realistic fluids with equations of state, such as p = ae where the 
pressure p is proportional to the energy density e, several important results have 
been obtained so far. There are self-similar collapse solutions which result in naked 
singularity formation |H]. The generic appearance of naked singularities has been 
strongly suggested numerically for a very small value of a [01El|. The critical collapse 
appears at the black hole threshold for < a < 1, which can be identified with naked 
singularity [TT1 IT2] . In spite of these findings, a comprehensive understanding of the 
endstates of gravitational collapse in this system has not been achieved yet. This is 
partly due to the fact that we do not know a general exact solution for this system, 
unlike for the dust cases. It should be noted that there are several related papers to 
this work by Christodoulou [13J. 

Recently, Giambo et al. gave an interesting approach to the spherically symmetric 
system of a perfect fluid with the barotropic equation of state p = p(e) [HJ HHj. They 
adopted the so-called area-radial coordinates and derived a second-order quasi-linear 
partial differential equation (PDE) for a mass function, which is identical to the Misner- 
Sharp quasi-local mass. Due to the presence of pressure the mass function is time 
dependent. A solution to this PDE gives a solution to the whole set of Einstein's 
field equations. Using mass functions which are analytical around the centre, they 
have examined the spacetime singularities resulting from the gravitational collapse. In 
particular, they claimed that there are always naked singularities for that class of mass 
functions. 

In this paper we analyse the behaviour of the metric around the regular centre and 
show that analytic mass functions imply staticity around the centre. This means that, 
in this case, a central singularity can not be formed without an infinite discontinuity in 
the fluid density. We then present a more general analysis of the final fate of collapse in 
this system by (i) studying other classes of mass functions and (ii) investigating naked 
singularities formation within such classes. 

The plan of this paper is as follows. In sectional we briefly recall the formulation of 
spherically symmetric spacetimes in area-radial coordinates. In section 0] we clarify the 
behaviour of metric functions during the regular evolution, which is a key to the present 
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problem. We investigate what type of collapse can be represented by the analytic mass 
functions given by Giambo et al. [T^ ITo] in section 0J In section we investigate the 
Friedmann solution and construct other dynamical mass functions. We investigate the 
relationship between the mass functions and naked singularity formation in section |3 
In sectional we give concluding remarks. We adopt the units in which G — c — 1. 



2. Formulation 

2.1. Comoving coordinates 

The line element of a spherically symmetric spacetime in comoving coordinates can be 
written as 

ds 2 = -e 2v dt 2 + e 2X dr 2 + R 2 (d6 2 + sin 2 Odcf) 2 ), (1) 

where v = u(r,t), A = A(r, t) and R = R(r,t) are of class C 2 . We consider a perfect 
fluid as a matter field and denote the energy density and the isotropic pressure of the 
fluid by e and p, respectively. Then, the Einstein field equations can be reduced to 

^' = 4ireR 2 R', (2) 

4/ = -4irpR 2 R, (3) 

R> = Rv' + R'X, (4) 

p' = -(e + P y, (5) 

where a dot and a prime denote derivatives with respect to t and r, respectively, and 
\J/ = \l/(r, t) is the Misner-Sharp quasi- local mass, 

*( r , t) = | [l - R' 2 e~ 2X + R 2 e~ 2v \ . (6) 
The matter (number) density is given by 

P= ^ER^ (7) 
where E is an arbitrary positive function of r. Here we assume a barotropic perfect 
fluid with a linear equation of state 

p = ae, (8) 

where a > is constant. This together with p = pde/dp — e implies e = p a+1 , where the 
constant of proportion was absorbed into the function E. Therefore from equation (0) 
we get 

e u = p~ a , (9) 

where the scaling of t is chosen so that the above equation holds. Using the freedom in 
scaling of r, we can set at the regular initial moment 

R(r,0) =r, (10) 

#(r,0) = l, (11) 

if there is no apparent horizon initially. In the following we choose this scaling. 
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2.2. Area-radial coordinates 



Using area-radial coordinates (r,R) we find from equations (j2j), (JHJ) and (JBJ) that both 
R' and the energy density are related to the mass function as follows 

a ^ r 



R'= - 



P 



a+l 



(12) 
(13) 



AiraR 2 ' 

where a subscript denotes the partial derivative with respect to the indicated variable 
in the area-radial coordinates (r, R) . In these coordinates, the line element (0) can be 
written as 

2 



ds z 



R' 



dR 2 - 2R'dRdr + — 1 - — ) dr 



Y 



2^\ 



R J 



where 



+ R 2 {dd 2 + sin 2 
Y = R'e- X , 



u 



VY 2 

R 



1. 



(14) 

(15) 
(16) 



Using equations (jJJ), (|T2^1 and (fTHJ) . the function K can also be written in terms of the 
mass function as 



Y = i : • ( 17 ) 

(a + l)p a v ; 

We can now use equations (0), (©, (|T3*j) and (JT7J) in equation (j^J) to reduce the Einstein's 
field equations ©-© to the following PDE 
1 C2$l 



a + l 



a)Y' 



^,RR _ ^,rR 



2^ 



a + 1 V R 



— -l + (l-a)F 2 --2r 2 - + — 



2a 



a+_l j^ y2 ^, rr 



a * 



, a 



R 



R R 2 aR 



0. 



a * >r £ 

Now, by eliminating E(r) from equation ()18|l. using two arbitrary functions ^o(r) 
ty(r,r) and Y Q (r) = Y(r,r), one can obtain the following PDE for 



2^ 



a + 1 V R 
2a /2tf 



- l + (l-a)F< 



a + 1 V R 



+ 2Y 2 - 1 



1 



I 1 — 

R R? aR 



+ 



a + 1 tt, fl 



1 ttg 



2a 



and 



F(r, fl) 



Mr, fl) 

^,r( r ) r ) 



a + 1 % (a + l)r 
^ jR (r,r)i? 2 l ° TT 



r 2 = o, 



a * 



Y n (r) 



(19) 



(20) 
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In the following, we call the above quasi-linear second-order PDE Giambo-Giannoni- 

Magli-Piccione (GGMP) equation. Because we adopt the scaling of the comoving 

coordinate r which is given by equations (fTT)|) and (JUJ, the mass function must also 

satisfy equation (|T2"|) at the initial surface, i.e. 

01 

*iir,r = — tf r (r,r). (21) 

a + 1 

3. Behaviour of the metric tensor around the regular centre 

First we consider the comoving coordinate system. The regularity condition at r = 
requires 

fl(0,t) = 0, (22) 

#(0,t) = e A(0 - t) , (23) 

during the evolution. It should be noted that these conditions are used in Giambo et 
al jUj. Also, in the following analysis, it is sufficient to assume lower differentiability, 
C 1 class, for the functions R(r,t), R'(r,t), e A(r,i ) and p(r,t) on i = const slice. 
An estimate of equation (0) at t — with a regular centre yields 

E(r) « / : „ , (24) 

V ; 4vrp(0,0)r 2 ' V ; 

where we have used (fTUj) . (|TT|) and (J23*)l and the weak equality is used with the following 
meaning: 

A(r)^B(r) & lim^ = l. (25) 
r^o B(r) 

Then, substituting equation (J2l*Jl into equation (JJJ), and taking condition (J23*|) into 

account, we get 

,2 

p(r,t)ttp(0,0)j^, (26) 

where the both sides depend on r and t and the weak equality is used with the 
following meaning: 

A(r,0«S(r,t) lim^4 = l, (27) 

r^o B{r, t) 

i.e. the limit is taken as t is fixed. Integrating equation for small r as t is fixed, we 
have 

a *"S^ +c w- (28) 

where C(t) is an arbitrary function. From the regularity condition (J22j) . we conclude 
that C(t) = 0. Then, we finally get 

R - "4541 r, (29) 




#rs cp^i . (30) 
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This implies that the area radius R is proportional to the comoving coordinate r in 
lowest order as t is fixed. It should be noted that the above discussion does not prove 
uniform convergence. 

Next we move on to the area-radial coordinate system. The above condition can 
be translated into the area-radial coordinates. A family of t = const regular spacelike 
hypersurfaces in the comoving coordinates is transformed to a one-parameter family of 
curves R = Rt(r) in the rR plane, which cross the origin R = r = along straight 
lines R = rr with < r < oo. If the fluid is static, the family of t = const spacelike 
hypersurfaces is transformed to one degenerate straight line R = r in the rR plane. If 
we consider models in which the central density changes monotonically with time, this 
slope parameter r parameterises the family of evolution curves. Moreover, if we focus on 
monotonically collapsing situations, we can assume < r < 1. The initial hypersurface 
t = corresponds to the straight line R = r with forty five degrees in the rR plane. As 
time proceeds, the slope of evolution curves at the origin R = r = gradually decreases 
towards zero. The r axis corresponds to the singularity curve. This means that we 
can approach a central singularity along a curve R ~ r n with n > 1. The regularity 
condition (J23|) is now translated into the area-radial coordinates in terms of the function 
Y as 

lim Y(r,rr) = 1, (31) 
for < t < oo, where the limit is taken as r is fixed. 



4. Analytic mass functions 

In this section we would like to investigate what type of collapse can be represented by 
the mass functions given by Giambo et al. [TH EH] . We recall that in jHj it is assumed 
that \I/ is an analytic function with respect to r and R which can be expanded as 

oo 

*(r,i*) = £ E *<^#, (32) 

k=0 i+j=3+2k 

where i and j are nonnegative integers, and it is proved that the regularity condition 
(13 lj) requires the following structure for the lowest order terms: 

*(^) = Kr 3 --^ 3 ) + E £ ^ riRi - (33) 

In addition, it is assumed that the coefficient ^41 vanishes in the proof of Lemma 3.5 
of Using this form of the mass function, it is claimed that the central singularity 
is naked in Theorem 3.4 of |14j . In the following, however, we show that their solution 
has a static centre in the regular evolution. 

Let us consider that we have regular evolution after the initial moment. During the 
regular evolution, we can use the relation R rs rr (0 < r < 00) for small r, as discussed 
in section 01 Then, it can be shown that the central density is 

limp (1+fl) (r,rr)= 3 \ (34) 
r-*o Snia + 1) 
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by substituting equation into equation (|T3j). that is, the central density does not 
depend on r. This guarantees, through equation 0, that the function e u is constant 
with time at the centre. This also guarantees, through equations ()29|) and (J3Uj) . that 
R k, r and i?' ~ 1. We can also obtain the same result from equation (|T2|) . Equation 
fT2|) with the relation R = rr (0 < r < oo) implies 



(35) 

Since this relation is satisfied at any time t, this can be integrated as t is fixed. The 
result is 

R 3 (r,t) « r 3 + L>(t), (36) 

where .D is an arbitrary function of t. The regularity condition (J22|) requires -D(t) = 
and therefore i? ~ r and i?/ ~ 1. Hence the regularity condition (J23*j) implies 
e A ~ 1. Therefore, the symmetric centre can not become singular without an infinite 
discontinuity in p. 

To construct mass functions which can describe genuine dynamical situations, we 
shall include terms of negative power of R in an expansion of the mass function. In the 
next section we shall construct such mass functions. 



5. Mass functions involving fluid dynamics 

We have pointed out some problems associated to the choice of mass functions (}3*3*j) . 
In this section, we construct other forms for \1/ which can be used to study models of 
gravitational collapse as well as cosmological models. 

To begin with, we consider the Friedmann solution as the simplest solution of a 
dynamical fluid. In this case, the matter density is homogeneous and satisfies 

poca(t)- 3 = 3 , (37) 

where a is a scale factor. Substituting this relation into equation (JT3J) , we obtain, 



Po - 

r 



R\~ 3 ( V 



1 

r J , 3f 



AnaR 2 ) 

where po is a constant. After integrating equation with respect to R, the mass 
function can be written as 

*(r,R) = ^ ^sr + ZCr), (39) 

where f(r) is an arbitrary function. From the assumption that R = r is the initial 
regular slice, we get that R = rr is also a regular slice, where < r < 1 is a constant. 
When we adopt the regularity condition for we find /(0) = 0. Furthermore, from 
equation (J2Tj) . we obtain /' = and then 

f(r) = 0. (40) 
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Therefore, the mass function (J39|) of the Friedmann solution results in 

r 3(a+l) 

*(r> R ) = *3-^T> ( 41 ) 

where # 3 = 4vrp^ +1 /3. 

In order to recover the familiar form of the Friedmann solution, we note that the 
metric functions in comoving coordinates e v and e A can be written in terms of the scale 
factor a using equation (fTHjl . Y = Y (r) and e u oc p~ a . Then, after an appropriate 
coordinate transformation we get 

ds 2 = -dt 2 + a 2 (t) \4tf^ + r2 ( de2 + sin2 ed( t> 2 ) J • ( 42 ) 

Next, we would like to investigate more general situations. A physically reasonable 
mass function should satisfy not only the PDE (j!9|) but also relation (|2*Tj) on the initial 
slice as well as the regularity conditions at the regular centre on the regular slice R = rr. 
In addition to these conditions, we can impose condition (A) that the mass function is 
the same as in the Friedmann solution (|4*T]) at the lowest order around the symmetric 
centre. To satisfy condition (A), the generalised mass function is written as 

y(r,R) = y/-^ + Il(r,R), (43) 

where Il(r, R) is a higher-order term on the regular slice, related to the inhomogeneity 
of mass distribution. The functional form of Il(r, R) is important for the final fate of 
the collapse, that is, whether it is naked or covered. 

In addition to the above restrictions on the mass function, we may impose condition 
(B) that the mass function recovers an appropriate dust limit as a — > 0, in which the 
mass function depends only on the comoving radial coordinate r. We can consider that 
the mass function \l/ is given by the following series: 

* = E^+^ 2fe+1 (-) > ( 44 ) 
k=i ^ r ' 

where {^2k+i',k = 1,2, •••} are constant coefficients and {fck+i'jk = 1,2, •••} are 
constants which depend only on k and a. Equation (|2*Tj) yields fok+i = —(2k + l)a. 
This satisfies both conditions (A) and (B). From this form of expansion, we can conclude 
that \I/ = ^ (r a+1 / R a ) . This together with equation (jl2j) yields R = r(t)r, where r is an 
arbitrary function of t. This implies that the singularity is simultaneous and therefore 
covered. 

A potentially more interesting form for the mass function, results from a slight 
generalisation of 



* = EEW,r 2fc+1 (^) , (45) 

k=l i=l v ' 7 

where ^ = {^2k+i,i]k = 1,2, = 1,2, •■■} is a constant coefficient matrix and 
{/32k+i,i',k = l,2,---,i = 1,2, •••} are constants which depend only on a, k and i. 
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Condition (A) yields /3 3 i = —3a and /3 3i = for i > 2. Condition (B) is satisfied if and 
only if fak+ij "~ > as a — ► 0. Equation (}2"T]) yields 

oo 

E *2*+i,i l#»*+i,t + ( 2k + !)«] = °- ( 46 ) 
i=i 

This implies that if the coefficient matrix \I> is regular, 02k+i,i — — (2A; + l)a is obtained, 
which eventually results in the same situation as the simplest model. Therefore, if the 
mass function which admits the expansion ()45|) describes naked singularity formation, 
then the coefficient matrix St must be singular, or det \I/ = 0. In this case, there exists 
at least a set of values for k and i such that fok+ij 7^ —(2k + l)a. 

We note that while the emphasis here has been given to models of collapse, the 
mass functions ^ can also be used in cosmological models. In particular, mass functions 
(jH| and (jH| can be associated to a simultaneous big bang and can be used in the study 
of the early universe. 



6. Naked singular mass functions 

In this section we investigate mass functions which involve spacetimes with naked 
singularities. We follow a technique used by Mena and Nolan |5, and developed by 
Giambo et al. JE]- To show the existence of null geodesies emanating from the central 
singularity, we show the existence of the region from which future-directed outgoing 
null geodesies cannot get out when they are traced back to the centre. This region is 
between two curves R = Ri(r) and R = R2^r) (0 < R±(r) < i?2(^)) both emanating 
from the central singularity. Defining (p(r,R) as the first-order derivative dR/dr along 
a future-directed outgoing radial null geodesic, i.e. 

^R) = -^^(l-l), (47) 



the curve Ri satisfies, 
dR x (r) 



dr 

while i?2 satisfies 
dR 2 



>f(r,Ri). (4* 



<ip(r,R 2 ). (49) 



dr 

In the following analysis, we assume that only one term in the expansion of Il(r, R) in 
S3), say 

n r 5 - m iT\ (50) 

where m is a constant, is important to examine the existence of naked singularities. 
Under this assumption we will derive the conditions on Ho and m for the existence 
of radial null geodesies emanating from the singularity which forms from gravitational 
collapse. As we have already mentioned, to construct a solution we may need another 
term 

U ir 5 ~ n R n , (51) 
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which is of the same order as the term of equation (|5(J|1 on the regular slice. 

In this analysis, an important role is played by the apparent horizon R = Rh{r), 
where 

R h (r) = 2*(r,R h (r)). (52) 

Under the above assumption, R h behaves around the centre as 

Rh = hr n \ (53) 

where the index rih > 1 depends on the value of m and h is a positive constant. After 
a careful analysis, we obtain rih and h in terms of the value of m, 

n h = 3 ( a + 1 \ h= (2^3)3^1 for m>-6a-l, (54) 

3a + 1 

3(a + 1) / / — n \ 

n h = 3a + 1 , h = (#3 + V^l + 2n J for m=-6a- 1, (55) 

n h = , h = (2I1 ) T ^ for m < —6a — 1. (56) 



1 — m 

It can easily be shown that 



d R 

v(r,R h ) = 0<—± = n h hr n »- 1 . (57) 
dr 

Therefore we can identify Rh with R^. 

If there is a curve R = R x (r) such that 

dR 

(a) R x (r)>R h (r) and (b) — - < <p(r, R x (r)), (58) 

dr 

we can identify R x as R2 and show the existence of radial null geodesies emanating from 
the central singularity. Here we assume that the curve R = R x (r) behaves as 

R x = xr n * (59) 

with n x > 1, which assures that R x emerges out of the central singularity. Condition 
(a) in equation (|58|l can be rewritten as 

(al) n x < rih or (a2) n x = rih and x > h. (60) 

In the rest of this section we would like to investigate the behaviour of R'(r, R x (r)) along 
the line R x , 

_ 3a(a + 1)* 3 + a(5 - m )Yi QX 3a+rn r^ a+m >-- 3a - m+2 
^' ~ 3a(a + 1)^3 - (a + i) m U x 3a+m r( 3a+m ) n -- 3a - m+2 

x xr^" 1 (61) 

and the limit, 

to ( 1 " = C (62) 



where C x is a constant for each x. Then, tp(r, R x (r)) can be represented as 

V (r,R x (r))^C x R'(r,R x (r)) (63) 



Spherically symmetric perfect fluid in area-radial coordinates 



11 



around r = 0. 

When the inequality 

3a + m — 2 < (3a + m)n x (64) 

holds, the first terms of the denominator and the numerator are leading terms and then 
the right-hand side of equation (JBTj) is approximated as R'(r,R x ) « xr nx ~ x . Also, it 
can be shown that Y(r, R x (r)) ~ 1 around the centre. Therefore, outside the apparent 
horizon, it can be easily shown that the limit C x is less than or equals 1. Therefore the 
following relation holds 

<p(r, R x (r)) ss C x xr n *~ x < n x xr n *~ x = (65) 

As a result, condition (b) in equation ()58|) cannot be satisfied and the central singularity 
would not be naked in this case. 
Next we consider the case 

3a + m — 2 = (3a + m)n x , (66) 

where R'(r, R x ) behaves as 

lr ' x> ~ 3a(a + 1)^3 - (a + l)mU x 3a + mXr ' 1 ' 

around the centre. Also, the following relations hold around the centre, 

^ _|_ (5-m)n x 3a+ m 

Y(r,RM)~- n ° +im ; zsz (68) 







"trip 3a+m ] "+ 1 



W(r,R x (r)) ^ 2n r 3(a + l)-(3a+lK 



- , .. ,. -Tox^ 1 r^-^ +i ^. (69) 



When m > —6a — 1, condition (a) means 



n. < 1 + (70) 



Using this inequality and equation (j66|) . it can be shown that 

2 

m = —3a < —6a — 1, (71) 

n x - 1 

which is inconsistent with m > —6a — 1. 

When m = —6a — 1, we obtain from equation (|Blj|) 

3(a + l) , , 

n x = -z — —j- = n h . 72 
3a + 1 

To satisfy condition (a2) the following relation should hold 



x > h = ( * 3 + V^i + 2n ) 3Q+1 . (73) 
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3a + 1 



< n 



\ 



i 



1 - 2^ 3 x 



-3a-l 



2n x 



-6a-2> 



(6a+l)n 
3a*3 



2a 

-3a-l 1 a+1 



1 _|_ 2n o t-— 3a— 1 
1 ~r * 3 X 



X 



i i 2TI a „-3o-l 

1 + ga X 

1 I (6a+l)n „-3a-l 

3«'I'3 



(74) 



To satisfy this inequality, the right hand side of it should be larger than unity. Therefore 
we need 



1 _|_ 2rio — 3a— 1 

^3 



1 _|_ (^ a + 1 ) u o x -3a-l 



> 1 



at least. This inequality holds when Ho < and 
„ 3a+ i . (6a + l)n 



X 



> 



3a^, 



(75) 



(76) 



We note that the positive matter density condition ^/ ^ < and the non-shell crossing 
condition \I/ ir (r, R) > 0, which are derived from equations (fT3*|) and (jl2J) . are satisfied in 
these cases. Now we can rewrite equation (f?4"|) as 



< 



- (l - 2^ 3 x 

2 3(a + l) 
3a + 1 



-3a-l 



- 2U x 



-6a-2 



(6a + l)n 



1 + ^ x -3a-l 
*3 



3a^ 3 

2 



x -3a-l 



a + 1 



+ 



• 3(a + l) V ( l (6a + l)no ^ 3a _; 
3a +1 / 1 3a^ 3 



• (77) 



This inequality can hold, for example, if a < 1 and x 3a+1 is sufficiently near its lower 
limit — - g^ 11 " given by (fTo^) . To be compatible with equation (f?3^) . it may be needed 
an additional condition for Ho as 
6a(3a + 1) 



< Iln < 0. 



(78) 



(6a + l) 2 

When m < —6a — 1, it can be shown that condition (al) is satisfied but (a2) is not 
from equation (JoT)j) . In this case Condition (b) gives 



3(a + l) 
3a + 1 



< 1 



\ 



mU Q 3a+r 
3a* 3 



2a 
a + 1 



1 _i_ (5-m)n 3 a+m 
1 ^ 3(a+l)*s 



-, _|_ (5-m)n 3a+m 
1 ^ 3(a+l)^ 3 

1 3a* 3 



(79) 



Similarly to the previous case, we can show the existence of the parameters set which 
satisfies (J79|) when a < 1 and n < 0. In this case, n does not have a lower limit. 
The last remaining case to analyse is 

3a + m - 2 > (3a + m)n X} (80) 

where the second terms of the denominator and the numerator of equation (j61)l are the 
leading terms. Around the centre, R'(r,R x ) behaves as 



a 5 — m 
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Along the curve R, 



■XI 




r 



-3a+2-m+(3a+m)nx 

CK + l 



oo (82) 



(83) 



and we conclude C x = for the curve which satisfies the condition (a). This means that 
the condition (b) is inconsistent with (a). 

In summary, there are parameters sets of m, ^3 n , and x which satisfy m = — 6a— 1 
and equations (|75jl and (j7j| or m < —6a — 1 and equation ((79)1 . In these cases, there 
are curves R = Rh(r) and R = R x (r), where Rh satisfies equation (|57j) and R x satisfies 
equation (jHHj) . Therefore when the mass function with such parameters is a solution of 
the GGMP equation, for which there is a central singularity after a finite time, then we 
can show the existence of radial null geodesies which emanate from central singularity. 

We note that the present analysis does not fully prove that the spherical collapse 
of a perfect fluid (with any a) produces naked singularities since in order to do that 
one should check, for each case, the consistency between the expanded form of mass 
function and the GGMP equation. After that we may use the above results to obtain 
a restriction on a which produces a naked singularity. Another way to do this would 
be to investigate whether the solution of equation (fTT?|l admits an expansion where the 
second-order term is given by equation ()50|) with m < —6a — 1 and the existence x 
which satisfies equations (fTHj) and (|73j) or equation ((75)1 . However, in general these 
tasks are much more difficult. As a final remark, we also note that while using the 
above procedure it should be checked, for each whether the behaviours of Rh and 
<p(r,R x (r)) are, in fact, determined by the second-order term of the mass function. 

7. Concluding remarks 

In this article we have investigated the collapse of spherically symmetric perfect fluids 
with the equation of state p = ae using area-radial coordinates. In this coordinate 
system, the physical variables are written in terms of the mass function which obeys 
a second-order PDE, which we call the GGMP equation. We have found that the 
assumption of analyticity of the mass function has the problem of implying staticity 
around the regular centre. In this case, the only way to form a central singularity is to 
permit an infinite jump of a physical variable. Furthermore, in order to have genuine 
dynamics, we argued that the expansion of the mass function around the centre should 
have negative powers of R. 

We have then constructed other forms for the mass function, starting with the 
Friedmann solution as the simplest dynamical solution of the GGMP equation. We have 
then considered classes of mass functions which represent inhomogeneous dynamical 
fluids. Such mass functions approach the Friedmann solution at the lowest order around 
the symmetric centre. We have also considered a condition for which the solution has a 
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meaningful dust limit. We note that the power indices and the coefficients of the mass 
function expansions were not given a priori and should be determined, in each case, by 
the GGMP equation. 

Using generalised forms for the mass function, we have investigated the conditions 
for the appearance of naked singularities resulting from collapse. We have proved that 
if (i) \1/ is a solution of the GGMP equation and is given by (|43j) with the second-order 
term proportional to r 5 ~ m R m with m < —Qa — 1, for which there is a central singularity 
after a finite time and (ii) there are x which satisfy equations (|7H|) and (|74jl or equation 
(J75J) , then there are naked singularities. One should be cautioned, however, that this 
result does not fully prove the existence of naked singularities. The work towards a 
complete proof of this is in progress. 
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